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£\\ • Abstract 



We analyse an asymptotic low-dimensional model of anti-plane shear in a thin 
i-Q • bi-material strip containing a periodic array of interfacial cracks. Both ideal and 

non-ideal interfaces are considered. We find that the previously derived asymp- 
r^ , totic models display a degree of inaccuracy in predicting standing wave eigen- 

frequencies and suggest an improvement to the asymptotic model to address 
this discrepancy. Computations demonstrate that the correction to the stand- 
ing wave eigenfrequencies greatly improve the accuracy of the low-dimensional 
model. 



1. Introduction 



> 

00 

l/»v . In this paper we present a method to correct discrepancies which arise in the 

r • asymptotic approximation of standing wave eigenfrequencies in a thin waveguide 

£— n . containing cracks and different types of interface. 

Substantial interest in the analysis of waves interacting with waveguide 
boundaries can be found in the literature. In acoustics and water waves, prob- 
lems in periodic waveguides have been studied in [5|, I 111 . |9[, among others. Bi- 
• »"H . material structures are widely used across many engineering disciplines, ranging 

p% ' from film coatings to armour production. In such applications, it becomes vital 

to understand how waves propagate through such structures and to estimate the 
stresses near the crack tip which may be sufficient to cause defects to propagate, 



which in turn may cause failure of the entire structure 12[ . 

The effect that interfacial cracks and other defects have upon the behaviour 
of structures is a particularly active research area attracting significant atten- 
tion. Modelling of interfacial cracks was studied in the important early papers 
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|Zi ll2|- The study of different types of interfaces, including ideal and non-ideal, 
is also widely covered in the literature, for example in [S, |j, la, IM llfll > among 
others. Recently the interaction of an interfacial crack with small impurities 
has been considered in the asymptotic regime in [14 L 

The present paper builds upon the results of 13[ and [16| and makes a 
breakthrough in improving the results in a wide range of cases that are impor- 
tant for applications. The former manuscript uses a weight function approach 
to construct an asymptotic model for out-of-plane Bloch Floquet wave propa- 
gation in a thin bi-material strip with an array of cracks positioned along the 
join. The latter considers a similar geometry, but with non-ideal interfaces lying 
between the cracks. Both papers use a weight function approach to obtain con- 
stants describing stress distribution near the crack tips and also dervice junction 
conditions for an asymptotic low-dimensional model. The conditions enable dis- 
persion diagrams to be constructed. Comparison with finite element simulations 
demonstrate that in both the perfect and imperfect interface cases, the low di- 
mensional model has high accuracy in the cases of the waves that propagate 
through the strip, usually of the order of 10 _4 %, but display a larger discrep- 
ancy in the case of the standing waves. The size of this discrepancy depends 
greatly upon material and geometrical parameters but is typically somewhere 
in the region of 3-15%. It is this standing wave discrepancy which we aim to 
address in this paper, by considering the next asymptotic terms for the solution 
and eigenfrequencies. 

The structure of the present paper is as follows. In Section [2] we present the 
problem formulation before developing the low dimensional model in Section 
[3l The eigenvalue correction term is derived in Section |4l firstly for an ele- 
mentary symmetric and homogeneous special case which is easily traceable and 
enables us to derive the correction term for the standing wave analytically, be- 
fore considering the fully general case. We then present numerical results which 
demonstrate the effectiveness of the correction method in a variety of cases and 
comment on situations where the model fails to give useful information. 

2. Problem Formulation 

The full problem formulation is as given in [l6| . but is summarised here. The 
geometry considered is a bi-material strip composed of two materials of shear 
moduli fix and /i2, with respective thicknesses eH\ and ei?2- The elementary 
cell of the periodic structure is shown in Figure [1] Along the interface lies a 
periodic array of cracks of length Z, between which the interfaces are imperfect 
whose extent of imperfection is described by the parameter k. The problem is 
singularly perturbed and so the case k = which corresponds to the perfect 
interface case requires different analysis [13J to the imperfect case [161 ]. The 
distance between adjacent cracks is a — I. The functions u^ (x, y) : j = 1,2, are 
respectively defined above and below the interface as solutions of the Helmholtz 
equations 
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Figure 1: Geometry of the elementary cell. 

A zero stress component is imposed in the out-of-plane direction along the top 
and bottom of the strip, as well as along the face of the crack itself: 



agHx^H,) = 0, a$(x, -eH 2 ) =0, i£ (-a/2, a/2), 
a$(x,0+)=0, 4 2 ;(x,0-) = 0, xG (-1/2,1/2). 



(2) 
(3) 



Outside the crack, along the boundary between lie and H E , the interface is 
described by the condition 

u^(x,0 + )-u^ 2) (x,0-) = eko-$(x,0 + ), xe (-a/2,-l/2)U(l/2,a/2). (4) 

Note again that the case n = corresponds to a perfect interface, while n > 
represents imperfect interfaces. We also assume continuity of tractions across 
the interface 



a$(x,0+)=<j$(x,0-), x€(-a/2,-l/2)u(l/2,a/2). 



(5) 



Solutions sought, uV\ represent Bloch-Floquet waves, so at the ends of our 
elementary cell x = ±a/2 we have for j — 1,2 the Bloch-Floquet conditions 



u { ' j) {-a/2,y) 
*i j ](-a/2,y) 



e-* Ka u^(a/2,y), y € (-eH 2 , eH x ), (6) 

e~ iKa a^(a/2,y), y € (sH^sH,). (7) 



2.1. Asymptotic Ansatz 

In both [13] and [16], eigenfunctions u(x,y) are approximated in the form 



u(x,y,s) = ]>> fc ^ m (^)+£V«(z,y)) 

fc=0 lm=l 

+ (w A k) (X A ,Y) + W B k) (X B ,Y^}+R N (x,y,s), (8) 

with scaled co-ordinates Xa, Xb and Y introduced in the vicinity of the left 
and right vertices of the crack defined as 



X A = 



X — XA 



X R = 



x — xb 



Y= V ~. 



(9) 



(k) 

Here, v m represent solutions of lower-dimensional problems within limit sets 
Qq , j = 1,2,3,4. Xm — Xm(x,y,e) are cut-off functions defined so that 
Xm(x,y;e) = 1 in £l s and decay rapidly to zero outside Q s . They van- 
ish near the so-called junction points A and B (the vertices of the crack). The 
terms Wj^ and W B represent the boundary layers near A and B, and V m 
is the 'fast' change of eigenfunctions in the transverse direction in the domain 
fie ■ Rn is the remainder term in the asymptotic approximation. 

The asymptotic approach described is employed in both [13j and [16[ . When 

(k) 

solving the low dimensional model to find the functions v m , computations dis- 
play a discrepancy for standing wave eigenfrequencies between the low order 
model and finite element simulations of the problem; these computations are 
presented in Section [5] In order to address the discrepancies that arise in this 
asymptotic model, in this paper we further consider the square of the frequency, 
to 2 , as an asymptotic quantity, writing 

N 

- 2 = E^- ( 10 ) 

fc=0 

It is not immediately apparent a priori that this amendment will lead to a 
large correction in the approximations of eigenfrequencies, but we will later see 
that this allows us to solve the first order low-dimensional model which causes 
a significant improvement in the accuracy of the model in those cases where 
the zero order model displays large discrepancies. Interestingly, in cases where 
the zero order model gives high accuracy, the corrections are very small. As 
an example, in one case we consider in Section [5j the first-order correction 
method alters the frequency of the first standing wave (for which the zero-order 
model gives a significant discrepancy) by 11% of its zero-order value, while the 
propagating waves (for which the zero-order model displays high accuracy) are 
only corrected by 10 -4 %. 

3. Solution of low-dimensional model equations 

Since the boundary layers Wa and Wb (see ([8])) decay exponentially, we 
have that far from the crack tip, Xm = 1, 



N 
U ~ 

k=0 m=l 






Substitution of this expression into the Helmholtz equation (flj and comparing 
coefficients of terms in e k , k = 0, 1 respectively yields the two equations 

(^ 0) )" + ^- + ^ 0) =0, m = 1,2,3,4; j = 1,2, (12) 

(^r + ^ + 4^ + 4^=0, m = l,2,3,4; J = l,2. (13) 



Expression (TT21 corresponding to terms in k = is the same as before when ui 
was not treated as an asymptotic series, but (fT3f is new. Above and below the 
crack, that is for m = 2, 3, we therefore have that 



2T/(1) 



d 2 v r 



8Y 2 






m = 2,3, 



(14) 



where d 2 — c\ and ^3 = c 2l which after integration and application of the 



boundary condition 



dY 



r=o± 

2 



yields the equation 



u-', 



(^rW + £?-t&>(s) + Jv^\x) = 0, m= 2,3. 



(15) 



For ?n = 1,4 (where no crack is present), rearranging and integrating (|13p and 
applying the condition for continuity of tractions across the imperfect interface 
yields the equation 



(^ 1) )"W + ^ ) W + ^ 0) W = o, 



where 



d\ = cic 2 i 



/ H1H1+ H2H2 
HiHic% + ii 2 H 2 c\ ' 



For the zero order approximation, 

,2 



{v<V)'\x) + %-vV>(x) = 0, m= 1,2,3,4. 



(16) 



(17) 



(18) 



3.1. Junction conditions and crack tip asymptotics 

The asymptotic representation of u does not affect junction conditions on 
the first two levels of the approximation. Junction conditions for the zero order 
approximation have been derived in [16| and read 



vf' (x A ) = 4° (xa ) = vf \x A ) ; 4° ( x b ) = 4° ( X B ) = 4° (xb ) , 
along with the conditions for flux 

^H 1 {vf ) )'{x A )+^H 2 {vf ) )'{x A ) = {vnHi +^H 2 )(v{ 0) Y(x A ). 
^H l {vf ) )'(x B )+p L2 H 2 {vf ) )'{x B ) = (pim+^H^v^YixB). 



(19) 



(20) 
(21) 



Junction conditions for the first order approximation at the right hand crack 
tip are given for m = 2, 3, by 



W/ 



v£>(x b )=v?(xb) + (-1) 



\m+l 



M2#2 



fiiHi + ^ 2 H 2 



a N A{(v^y}(x B ), (22) 



where N = P if n = (the perfect interface case) and N = I if n > (the 
imperfect interface case), and 

A{(v^)'}(x) = (vi 0) )'(x) - (vi 0) )'(x). (23) 

We will continue to use this notation throughout the rest of the manuscript. 
The junction conditions (|2"2")l are valid for both perfect and imperfect cases, 
but the form of the corresponding constants ap and ai are absolutely different 
and come from different analysis; this arises from the fact that the problem is 
singularly perturbated and so different analysis is needed in the cases k > 
to the case k = 0. The definitions of the constants oti and ap are derived in 



and [13J respectively and stated here. For the perfect interface case, the 




constant is defined as 



ap = ^ im„ I f(t)dt — In 



where 

_ if* - ta,nh(tH*) coth(t) _ /xi - fi 2 „ _ H ^ ~ H ^ ( 2 k) 

/U ~ (sinh(i) + M * siiih(tF*))*' '"* ~ Mi + M2 ' * ~ Hi + H 2 ' l j 
In the imperfect interface case, the constant is given by 

CT , = (ffl + ff2) {l/!^ ( + J_}, (26) 

where 

t ( 2 t(l + H m ) 2 , *(1-F») 

9{t) = AfT^ ^ + K,(l + /i,) COth ^^ + K,(l - /X,) COth —T~ 

(27) 

x m , „ n l^lH 1 +fi 2 H2 K(Ml+M2) , oc s 

K = {H l +H 2 )\ — — — , tt*= ■ (28 

y [ixix 2 B x B 2 K, (Hi + H 2 ) 

We stress that ai is a constant that depends heavily upon k and so describes 
how the junction conditions are impacted by the imperfect interface. The first 
order fluxes satisfy the relationship 

(fiiHi + f M2H 2 )(v^ ) )'(xB) - pnHi(v^)'(x B ) ~ mHi^y&B) = 0. (29) 

The analogous conditions for m = 2,3, at the other vertex can be obtained 
by replacing B by A and m + 1 by m in equation (]22[) . The other crack tip's 
condition for fluxes is as in equation (|29p. but again replacing S by A. 

The zero-order and first-order constants describing the singular behaviour 
of the full solution near the crack tips derived for the perfect and imperfect 
interface cases respectively in [13j and [16| are unaffected by the consideration 
of u> as an asymptotic series. However, if one continues to deeper levels of the 
asymptotics, the junction conditions of fifth order and higher would be affected 
by taking uj as an asymptotic series. 



3.2. Corrected low dimensional model 
3.2.1. Zero order low dimensional model 

Solutions of the zero order equation (fT5)) for m = 1,2,3,4 are of the form 



-L 0) (*) = A™ sin \2L X j + fl£> eos (jj2L* j , (30) 

The first order equation ([TH]) has solutions in the form 

v£\x)=AW B m(^x)+Bg>c OB (£-x) + ufF m (x), (31) 

\«m / \d, m J 

where 

'-w-ss;^"^)-^ -(£•)}■ <32 > 

We note that assuming the zero order system has been solved, all constants in 
this expression for F m (x) are considered known. 

Let us first consider the zero order case. We see from (|3"0|) that eight con- 
stants need to be evaluated, A m and B m for m = 1, 2, 3, 4 which we write in 
the column vector A^ defining the notation 



A (k) 



A (k) B (k) A (k) B (k) A (k) B (k) A (k) B (k) 



(33) 



We have six junction conditions: two from ([T9|) and one from each of ([20]) and 
(|2~Tj) , which can be complemented by the two Bloch-Floquet conditions to yield 
the 8x8 matrix equation 

M 8xS A^ = 0. (34) 

We introduce notation to abbreviate the entries of Mg X 8 a s follows: 

S m = sin ( — i-rooXB I ; C m = sin ( — -w x B I ; m = 1, 2, 3, 4; (35) 

\«m / \d m J 

^i = — tj \ 3 u ~r~' 3 = 1 > 2 ; S a = sin ( w - ) ; C a = cos (^n-) , (36) 
fiiHi + n 2 H 2 dj+i V 2/ V 2/ 

where njj = Wj/d\ is introduced for normalisation. Now, M$ x % = 

-C 4 

-c A 






a „ — iKa/~i 

ar~i „ — iKa 









5' 2 


c 2 








— £4 














S3 


c 3 


— 04 








ikCt 


-i/'i'Sa 


1P2C3 


-^2^3 


-c 4 


-Si 


-Ci 


s 2 


C 2 











S'i 


-Ci 








-5 3 


c 3 





-Ci 


-Si 


V'iCs 


^1^2 


^2^3 


■02 5 3 





-5a 


c a 














p — iKa 


c a 


•J CI 














p —iKa 



(37) 



The determinant of Mg x s can be written in the form 

det(M 8x8 ) = A(Lo )e- 2tKa + B{u )e^ Ka + A{uq). (38) 

It can be shown that for the case in which all wavespeeds d m are equal to 
d, say, both A(uiq) and B(ljq) are zero when loq = mrd/(2xB), n £ N. It 
follows that in this case, the standing waves have no dependence upon the 
Bloch-Floquet parameter K. This observation motivates us to consider the 
special case discussed in subsection 14. II 

3.2.2. First order low dimensional model 

Applying the junction and Bloch-Floquet conditions for the first order equa- 
tion (|18|) yields the matrix equation 

M 8x8 v4« - w\N &x8 A^ + B A A{(v^Y}{x A ) + B B A{(v^)'}(x B ). (39) 

Here, M 8x8 is the matrix defined in (|37|) . A^ and A^> are the coefficients de- 
fined in (|33[) . Since M 8x8 , is singular, this equation gives a solvability condition 
which will allow us to find the correction term, u>i- The matrix N 8x8 is defined 
as 

Q Q — d\X B C 2 d^X B S 2 Q Q 18^4 -X B Sj 

1d2tOQ 1d2tOQ 2tuq 2cuq 







U7n 

-dtxsCs dtxsSs x B Cj -x B Sj 

2tl3t^7Q 2(l3tUQ 2tuq 2tUQ 



7V 3 ,3 #3,4 #3,5 N 3fi N 3J 7V 3 , 8 

-x B C\ -x B Sx d 1 x B C 2 d 1 x B S 2 g 

2tt7o 2roo 2d2^ T o 2c^2tUQ 

-x B C\ -x b S\ q q d\x B C-$ d\x B Sz g g 

2tuo 2tt?o 2d3^7o 2rf3tno 

#6,1 #6,2 #6,3 #6,4 #6,5 #6,6 

aC a aS a a ZC a —aZS a 



AzOq 4tx7o 4zuq Atuq 

#8,1 #8,2 iV 8 , 7 N 8A 

where Z = e - lKa and for q e {3, 4, 5, 6, 7, 8}, 



(40) 



^ = 5[^7 XB U 2 J C L^J + l 2 J S '-V. 

Forge {1,2,3,4,5,6}, 

(-i)« +1 di //i + (_i)«\ , A -(-i) 9 



^ M = d La±1J x Hl — 2 — J c wv — 2 — J'WJ 

and the expressions of the eighth row are given by 

N 8A = ^S a --^-C a , N 8a = -^-S a -^C a , (43) 

4 ZVJc) ZWn 4 



N> 



8,7 



„-iKa 



7, Co - -rS a ) , Ng,8 

ZWq 4 



-'--"" {^ s - + \ c - 



The vectors Ba and Bb are given by 



Ba = a N 
Bb = ctN 







fJ. 2 H 2 



M2#2 



HlHi 



HiH 1 +fj l2 H 2 fJ.iH 1 +fj, 2 H 2 







HiHi 



HiH 1 +fj, 2 H 2 HiH 1 +p, 2 H 2 







(44) 

(45) 
(46) 



4. Derivation of first order correction term, Wi 

4-1- Homogeneous Symmetric Case 

In this section, we condsider the symmetric case in which Hi = Hi and [i\ = 
/i2. This simple case is instructive since the symmetry enables us to analytically 
determine the eigenfrequency of the first standing wave in an easily traceable 
process; we will later make indications on how the method for the general case 
relates to and differs from this procedure. Moreover, this eigenfrequency does 
not have any independence as is the case for inhomogeneous setups, which 
enables us to easily separate the first standing wave solution from the others. 
In the case of the standing wave, the beams above and below the crack vibrate 
while the others do not (that is, v[ (x) = v\' (x) = for k = 0, 1). 

The problem formulation for the symmetric case is as follows. For the zero- 
order approximation, solutions satisfy (IT51) with d m — d for all m = 1,2,3,4. 
To isolate the standing waves (whose frequencies we wish to impose a correction 
upon) we impose the condition 



,((') 



(0), 



v { ;"(x)=v)l"(x)=0. 
The junction conditions for the zero order approximation then simplify to 



(47) 



V X 2 '{XA) =«2 ( X B) 



„(°) 



0, v$"{x a )=v$"(xb) = 



(0), 



along with 



(vi 0) )'(xA) + (v^Y(xa) = 0, (v^Y(x B ) + (vi 0) Y(x B ) = 0. 



(48) 



(49) 



For m = 2,3, we have that the general solution of of the zero order LDM is of 
the form (J30J) with d m = d. Applying conditions (05])-(|31]) yields 



(o), x D /wo \ (0)/ s 

i> 2 (x) — Bcos I — x J , v^'(x) 



= —Bcos 



(t*) 



(50) 



The first order approximation equation is of the form (1151) with d r , 



,(°) 



Since vin are now known functions, the corresponding system consists of the 

(51) 



two ordinary differential equations 

i 

d 2 



(vi 1) nx) + ^\x) + ^Bcos(^x)=0, 



d 2 




Figure 2: Dispersion diagram for the perfect interface case with medium crack length I = 2m 
in an elementary cell with a = 6m. The material above and below the crack is iron. The solid 
black lines show finite element results, while red crosses (x) show the zero order approximation 
and blue circles (o) show the corrected regime, with the first standing wave corrected by the 
analytic derivation of wi as presented in equation Il58|l . 



These ODEs have respective elementary solutions 

v£\x) = 4 1} sin (^x) + B« cos (^x) - uf 

vP{x) = 4 X) sin (^x) + B« cos (^x) + u\ 
complemented by the juction conditions (which follow from |16| ) 



0. 



Bxsm(^fx) 
2ujqcI 

Bxsm(^fx) 
2(jJq<1 



vSiHxp) 



(-1) 



m+p 



a N A{(v^y}{x fj ), m = 2,3; = A,B, 



where p = if /3 = A and p — 1 if f3 — B. For the first standing wave, w 
and so applying (|55|) with j3 — B, m = 2, we see that 

,(i) B ( 2 12 
Applying the condition (155)) with ft = A, m = 2, yields the condition 



A 



(i) 



d2 



' 47T 



- woa^v 



(52) 

(53) 
(54) 

(55) 
nd/l, 

(56) 
(57) 
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Since A 2 is equal to the positive and negative of the same function, it must be 
zero, from which we deduce that 

«l = p . ( 58 ) 

Thus, for the symmetrical case where both materials are the same, we have 
found an expression for the correction term <J\. 

Figure [5] shows the dispersion diagram for this homogeneous, symmetric 
case. The red crosses on the diagram indicate the zero-order approximation 
(u>o) of the eigenfrequencies, while the blue circles show the corrected first-order 
approximation. The black lines result from a finite clement computation. The 
derived correction method improves the standing wave frequency discrepancy 
from 3.7% to just 0.26%. 

4-2. General case 

We now consider the general case in which the materials above and below 
the crack and interface may have different thicknesses and shear moduli. For the 
first order approximation, after the application of junction and Bloch-Floquet 
conditions we obtain a matrix equation of the form 

MAW = NA^ + B A A{(v^)'}(x A ) + B B A{(v^)'}{x B ). (59) 



Here, M and N are both 8x8 matrices as defined earlier in (f37|) and (|40|) 
respectively whose elements depend on the Bloch-Floquet parameter K and the 
eigenfrequency ujq which is such that det(M) = 0. Because det(M) = 0, M has 
zero among its eigenvalues. We can write 

M = VDV- 1 (60) 

where V is a matrix whose columns are eigenvectors of M and D is a diagonal 
matrix with the respective eigenvalues of M along the diagonal. Premultiplying 
(|59|) by V' 1 , we can write 



V~ 1 MVV- 1 A < - 1) = V^iwlNA^ + B A A{(v {0) y}(x A ) + S B A{(i;( ))'}(^)), 

(61) 
which upon substitution of (|60j) becomes 

DV- X A^ = ufV- l NA^ 0) + V- 1 BxA{(v (0) ) / }(xa) + V- 1 B fl A{(t; (0) ) / }(a!B). 

(62) 
Since M is singular, it posesses zero as an eigenvalue, and so one row of the left 
hand side matrix in (1621) is equal to zero. Let us denote that row I. Then 

uj\ (y- 1 NA^) i + (V- 1 B A ) l A{(v^)'}(x A ) + (V- 1 B B ) l A{(v^)'}(x B ) = 0. 

(63) 
All matrices and parameters in this equation are now known, with the exception 
of cof which can now be written in terms of known quantities: 

(V-HB a A{(vW)>}(xa) + B B A{(v( ))'}(x B ))) l 

^ = (y-iAM(o)), ' (64) 

ii 



A potential problem with this computational method is that the matrix V 
may have a determinant which is close to zero when eigenvalues of M are close 
together. To eliminate any possible errors arising from this, we introduce a 
second computational scheme for computing w\. 

The Schur decomposition states that if A is a n x n square matrix with 
complex entries, then A can be expressed in the form A — QUQ^ 1 where Q is 
unitary and U is upper triangular, with the eigenvalues of A. In our case, we 
apply Schur decomposition to the transpose of M: 

M T = QUQ- 1 . (65) 

Schur decomposition is not unique; we may place the smallest eigenvalue in the 
first position along the leading diagonal, and since M is singular, this eigenvalue 
is zero. The first column of the upper triangular matrix U is therefore a row 
of zeros. Since M = (Q~ 1 ) T U T Q T ', we can premultiply the first order matrix 
equation by Q T and substutite to obtain 

V T Q T AW = u,lQ T NAW + Q T B A A{(v^)'}(x A ) + Q T B B A{(v^)'}(x B ). 

(66) 
The first row of the left hand side is zero, whence 

(Q T (B A A{(v^y}(x A ) + B B A{(v(Vy}( XB ))) 1 

l -. (67) 



(Q T NA«») 



5. Numerical results 



5.1. Materials and geometries used in numerical simulations 

For our numerical calculations we will consider a strip whose elementary cell 
is of length a = 6m with an overall thickness of H1 + H2 = 0.15m. This geometry 
corresponds to a value of e = 0.025. We will compare results from the low 
dimensional model against those from finite element simulations (COMSOL). 
We stress that finite element simulatons are efficient for comparison only in cases 
when the strip is not too thin, i.e. when e is not too small. The low dimensional 
model, however, remains valid as e — > 0. For our computations we vary four 
parameters as listed in Appendix |Appcndix B| These parameters are type of 
interface (perfect, imperfect, highly imperfect), length of crack (short, medium, 
long), materials (iron/aluminium [similar wavespeeds], magnesium/aluminium 
[less similar wavespeeds]) and thicknesses of each material (symmetric geometry, 
asymmetric geometry). 

We present in this section a number of dispersion diagrams, plotting fre- 
quency u> against the Bloch-Floquet parameter K. We refer to plots of u = loq 
as the zero order approximation, and to plots of u calculated according to (flQ| 
as the first order approximation, or the corrected frequency. 
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Figure 3: Dispersion diagram for a perfect-interface strip composed of equal thicknesses of 
aluminium and magnesium with a medium length crack (I = 2m). The solid black lines show 
the finite element results, while red crosses (x) show the zero order approximation and blue 
circles (o) show the corrected first order approximation. 



5.2. Perfect interface 

5.2.1. Materials with similar wave speeds 

Figure [3] demonstrates the effectiveness of the method of eigenfrequency 
correction for the standing waves. This dispersion diagram results from com- 
putations corresponding to the case of a strip with a 'sensible' crack length, 
composed of equal thicknesses of magnesium and aluminium (materials pos- 
sessing similar wave speeds). The diagram demonstrates that the zero order 
approximation agrees to a very high degree of accuracy with the finite element 
results in the cases of the waves which propagate through the strip (the slanted 
lines). Typically the zero order approximation for these waves' eigcnfrequencics 
differs from the finite element simulation only by around 10~ 4 %. However, it 
is clear that there is a discrepancy between the zero order model and the finite 
element results in the case of the standing waves (horizontal lines on the disper- 
sion diagram). The corrected first-order model retains the excellent accuracy for 
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Figure 4: Dispersion diagrams for a perfect-interface strip composed of equal thicknesses of 
aluminium and magnesium Left: Short crack (I = 0.6m) Right: Long crack (Z = 5.4m). The 
solid black lines show the finite element results, while red crosses (x) show the zero order 
approximation and blue circles (o) show the corrected regime. 



propagating waves, slightly increasing the accuracy while remaining of the order 
of 10~ 4 %, and hugely improves the discrepancy of the standing wave frequen- 
cies. While the correction is not completely uniform since the standing waves' 
frequencies depend upon K except in the case where materials have identical 
wavespeeds, a typical discrepancy for the first standing wave has decreased from 
3.9% to 0.38%. This can be considered as a surprisingly large correction, since 
there is no reason a priori to suspect that considering u> as an asymptotic quan- 
tity should cause such an improvement in the accuracy of the low dimensional 
model's approximation of the standing wave frequency. 

The dispersion diagrams for the cases with short and long crack lengths, 
again with materials of similar wave speeds, perfect interfaces and the same 
thicknesses of both material are contained within Figure^] As one would expect, 
the length of the crack does not significantly alter the eigenfrequencies of those 
waves that propagate through the strip (again the correction is on the level of 
10 _4 %), since they are not strongly influenced by the presence of the crack. In 
similar agreement with physical intuition, the first standing wave for the long 
crack is of much lower frequency than in the geometry housing the particularly 
short crack. The correction offered by the first order approach is relatively small 
in the long crack case, but since the zero order model in this case already gave 
good accuracy with only a 1.3% discrepancy for the first standing wave, this 
is not surprising. The corrected eigenfrequcncy of this wave agrees with finite 
element results to within 10~ 3 %. This can be seen as a surprisingly effective 



14 



correction since in the long crack geometry, the crack tips are close to the ends 
of the elementary cell. This gives the boundary layers surrounding the crack tips 
a small area in which to decay so that they do not influence the Bloch-Floquet 
conditions. 

In the case of the short crack, the zero order approximation of the first 
standing wave eigenfrequency is easily seen to be significantly different to the 
true value found in the finite element simulation, with a 12.2% discrepancy. 
After applying the correction method, the discrepancy decreases to 0.95%. 

5.2.2. Materials with more contrasting wave speeds 

The standing wave dispersion diagrams for a strip of aluminium and mag- 
nesium are presented in Figure [31 for both the symmetrical and asymmetrical 
cases. The correction is largest, as is the zero order discrepancy, in the sym- 
metrical case. The correction in the case of materials with different wave speeds 
is less uniform than in the Al-Mg case; this can be readily seen in the highest 
frequency standing wave shown in the left hand subfigure of Figure [5J The cor- 
rection still offers a significant improvement in most cases, however, although it 
is harder to quantify the exact size of a typical discrepancy. 





Figure 5: Dispersion diagrams for a perfect-interface strip composed of aluminium and iron 
containing a medium-length crack (I = 2m). Left: Equal thicknesses of iron and aluminium. 
Right: Different thicknesses (eHi = 0.01m thickness of aluminium, ei?2 = 0.14m thickness 
of iron). The solid black lines show the finite element results, while red crosses (x) show the 
zero order approximation and blue circles (o) show the corrected regime. 
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5.3. Imperfect interface 

The results in the case of the imperfect interface analogue follow broadly the 
same qualitative pattern as in the perfect interface case. The left hand subfigurc 
of Figure [6] gives the dispersion diagram for an iron- aluminium strip, joined with 
a thin layer of epoxy resin adhesive. Due to the different wavespeeds, the sizes 
of the standing wave corrections are dependent on the Bloch-Floquet parameter 
K, but in most cases the correction gives a significant improvement in accuracy. 
An interesting phenomenon can be observed when eigenvalues are close to each 
other in this subfigure; a zoomed section of the dispersion diagram to illustrate 
this is given in Figure [7] In Figure [3 some of the circles have been replaced 
by squares; these are the points which approximated propagating waves in the 
zero-order model which are corrected to approximate the standing waves for 
some values of K . In doing so, a crossing-over phenomenon occurs, where the 
order of eigenfrequencies switches after correction. The phenomenon becomes 
more pronounced at higher frequencies. 

The right hand figure of Figure [5] corresponds to a case where the materials 
are bonded in a highly imperfect fashion, using an adhesive with a shear modu- 
lus just a tenth that of the epoxy resin whose parameters are given in Appendix 
|Appendix A| (i.e. a material with shear modulus 2.5 x 10 8 Nm~ 2 and mass 
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Figure 6: Dispersion diagrams for imperfect-interface strips composed of equal thicknesses 
of aluminium and iron Left: Bonding material corresponding to epoxy resin (re* = 2.88). 
Right: A highly imperfect interface (re* = 28.8) representing an extremely soft bonding 
material. This rightmost subfigure demonstrates an example in which the low dimensional 
model is not effective. The solid black lines show the finite element results, while red crosses 
(x) show the zero order approximation and blue circles (o) show the corrected regime. 
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Figure 7: Zoomed section of the imperfect interface case, showing how the eigenfrequency 
correction method causes a crossing-over phenomenon. 



density 1850kgm~ ). This subfigure serves the purpose of presenting a case in 
which the low dimensional model ceases to provide entirely useful approxima- 
tions. Here the finite element simulation displays qualitatively different features 
which the low dimensional model does not predict at all. The fact that the low 
dimensional model breaks down is not surprising; in fact, the manuscript [iff 
derives a condition on k for the analysis on which the low dimensional model 
is based to be valid. The key point of this condition is that if k is too large, 
the boundary layers Wa and Wb (see equation © on paged) decay sufficiently 
slowly for the assumption that they are independent to cease to hold. More- 
over, if the boundary layers decay slowly from the crack tip, the Bloch-Floquet 
conditions will be influenced by the boundary layers. This case is interesting in 
itself and requires separate analysis. This phenomenon has been discussed in 
1- 



5.4- Conclusions 

The comparisons between the low dimensional model and the finite element 
simulations demonstrate that the proposed method of eigenfrequency correction 
is highly effective in most cases, typically improving accuracy for the standing 
wave eigenfrequency by an order of magnitude. Moreover, the correction it- 
self is computationally very efficient. The correction becomes less uniform for 
materials with significantly different wave speeds, improving the accuracy by 
different amounts in different parts of the dispersion diagram, and misses quali- 
tative features for setups with highly imperfect interfaces. Practically, however, 
such highly imperfect interfaces are unlikely to be encountered. Fracture pa- 
rameters are not affected by the analysis and for that reason are omitted in this 
paper and refer the reader to 13, [l6| , where discussions and details relating to 
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these parameters can be found. We only underline here that SIF (in the per- 
fect interface case) or COD (in the imperfect interface case) can be constructed 
as functionals on the low dimensional model without further work. For that 
reason, the results of this paper are of additional importance. 
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Appendix A. Material constants used for computations in Section [5] 



Material 


Shear modulus 


Mass dens 


ity 


Wave speed 




[Nm" 2 ] 


[kgm -3 


1 


[ms _1 ] 


Iron 


82 x 10 y 


7860 




3230 


Magnesium 


17 x 10 9 


1738 




3128 


Aluminium 


26 x 10° 


2700 




3103 


Epoxy resin 


2.5 x 10 9 


1850 




1162 



Appendix B. Details of parameters varied for computations in Sec- 
tion \E\ 

1. Type of interface 

• Perfect. 

• Imperfect — in the finite clement computations, a thin layer of 
epoxy resin is used. This corresponds to a value of k* = 2.88 in the 
asymptotic model. 

• Highly imperfect — in this case, the bonding material has a shear 
modulus a tenth that of epoxy resin. This corresponds to a value of 

k* = 28.8. 

2. Length of crack 

• Short — Crack length of / = a/10 = 0.6m. 

• Medium — Crack length of I — a/3 = 2m. This can be viewed as a 
'sensible' crack length. 
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• Long — Crack length of I = 9a/10 = 5.4m. 

3. Materials 

• Iron/ Aluminium — see Appendix |Appcndix A for shear moduli, 
densities and wave speeds. 

• Magnesium/ Aluminium — both materials have similar wave speeds. 

4. Thicknesses of each material 

• Symmetrical geometry — sHi = eH 2 = 0.075m. 

• Asymmetrical geometry — sH\ = 0.01m, eH 2 = 0.14m. 
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